APS 



Scaling of the Background Brain Dynamics and Alpha Ryhthm 

D.C. Lin 1 , A. Sharif 1 , H.C. Kwan 2 
1 Department of Mechanical and Industrial Engineering, 

IT) , 

O Ryerson University, Toronto, Ontario, Canada 

O ■ 

CN ; Department of Physiology, University of Toronto, Toronto, Ontario, Canada 

(Dated: February 2, 2008) 

Abstract 

The scaling property of the brain dynamics is studied based on the zero-crossing of the local 
electroencephalographic (EEG) recording taken from healthy young adults in eyes closed and eyes 
open. Evidence of coupling between the EEG fractal dynamics and the a rhythm is presented. 
An organization principle that governs this coupling relationship is proposed. In the a dominant 
q , brain state, a possible interpretation using the self-organized criticality similar to the punctuated 

equilibrium is discussed. 
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I. INTRODUCTION 



The cortical activity of the human brain in wakefulness and eyes-closed typically exhibits 
the 8~12 Hz a rhythm 1 . Although its origin remain open 2 ' 3 , this classical rhythm has been 
associated with the "resting" state of the cortex with most pronounced activity recorded in 
the parietal and occipital areas. 

While the a rhythm represents one of the major cortical activities, there typically lies 
a background "noise" component of unknown functionality and origin 4 . There is growing 
evidence of the importance of this background component. For example, local surface scalp 
measurements based on electroencephalography 4 ' 5,6 ' 7 (EEG) and magnetoencephalography 8 
reveal its fractal characteristics coexisting with moderate a rhythm. Similar results were 
found in multi-channel recording 4 ' 6,7 ' 8 , suggesting the global nature of the fractal dynamics 
in the cortex. The relevance of this background fluctuation in physiology may be seen from 
its state dependence property 8,22 where normal individuals in eyes closed average to a larger 
scaling exponent than in eyes open. Recent works further imply that the EEG fluctuation 
may have a broader implication on the other autonomic function such as the cardiovascular 
regulation 9,10 . 

The purpose of this study is to investigate the EEG background fluctuation and its 
coexistence with the a rhythm. While they are separately subjects of intense study, far less 
understood is the organization of these two prominent features of the brain dynamics. In 
addition, the findings reviewed above may require further clarifications when the a oscillation 
becomes dominant. For example, a trained meditator or Yoga practitioner can shift a 
significant portion of the EEG signal power to the a band that obscures any potential fractal 
characteristics in the background fluctuation. Fig. 1 shows EEG's with moderate and strong 
a rhythm from a normal subject and an experienced Yoga practitioner, respectively. It is 
observed that the power law trend diminishes almost entirely in the a dominant EEG, but 
not in the a moderate EEG. This apparently inverse relationship with the strength of a 
rhythm implies either the fractal background fluctuation is buried in the a oscillation, and 
can thus no longer be detected from the amplitude characteristics of the EEG, or there is 
simply no long-range correlated fluctuation in the a dominant brain state. Using the detrend 
fluctuation analysis 11 (DFA) on the integrated a dominant EEG actually supports the latter 
scenario of a white noise process in the background fluctuation (DFA exponent ~0.5); see 
Fig. lc. However, as we show below, it is not a valid description for the a dominant brain 
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state. 

In this work, we use the zero-crossing property to study the EEG background fluctu- 
ation coexisting with a rhythm. EEG zero-crossing has been primarily used to extract 
event-related frequency information 12 . Despite the obvious advantage of being insensitive 
to amplitude artefacts, the use of zero-rossing to study EEG background fluctuation ap- 
pears scarce. Watters and Martins used DFA to analyze the "EEG walk" constructed from 
the zero-crossing of EEG showing moderate a rhythm 13 . These authors found evidence of 
scaling and rejected the (uncorrelated) random-walk interpretation for the EEG background 
component. But we should point out that a DFA scaling exponent ~ 0.5 (random walk) is 
in fact possible in the a dominant brain state (Fig. lc). 

By excluding the a wave zero-crossing, we will show that the EEG background does ex- 
hibit fractal characteristics. For the scaling analysis in a dominant case, the zero-crossing 
approach is found to be more effective than other amplitude-based methods such as the power 
spectral density function and DFA. The complementary set of the background zero-crossing 
in the real line captures other EEG activities. For the a dominant case, this complementary 
set describes mainly the a dynamics. With the confirmed EEG fractal background, we con- 
jecture fractal scaling in the emergence of a dynamics. Support to this conjecture is found in 
the power law distribution of the a interval. Our analysis also suggests that the a dominant 
brain state may be interpreted in the unversality class of the self-organized criticality 15 of 
punctuated equilibrium 14 . While SOC has been proposed for the EEG fluctuation in the 
a frequency band 7,8 , and its analytical phase 4 , our result provides further evidence of its 
possible connection to the organization of the EEG background dynamics and the a rhythm. 
Finally, we are able to write down the governing principle for this organization that relates 
the EEG fractal scaling and the emergence of a dynamics. 

In the next section, the main idea of extracting the zero-crossing of EEG background 
fluctuation is described and verified with artificial examples. In section III, the method 
is applied to the EEG records taken from subjects showing little to significant a rhythm. 
Discussion and concluding remarks are given in the last section. 
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II. METHODS AND NUMERICAL EXAMPLES 
A. Main Ideas 

Let the EEG be x(t). The zero-crossing time is the level set {ti,x(ti) = 0} where the 
index % registers the order of the zero crossing event. In practice, is first determined 
by linear interpolation and then used to define the set of crossing-time-interval (CTI) C = 

= t i+ i — ti}. 

Zero-crossing of a stochastic process is a surprisingly hard problem; see, e.g., Ref. 16. 
Thanks to its self-similarity, the CTI for a fractal process is known to follow a power law 
distribution 17 : p(r) ~ r~ v , where p(r) is the probability density function (PDF). For exam- 
ple, v = 2 — h for the fractional Brownian motion (fBm) B^it) of Hurst exponent h. 

If fractal exists in the a dominant EEG, it can be captured in C\A a where A a denotes 
the CTI of the a oscillation. However, fractal CTI's can occasionally lie in the a rhythm 
range 18 . It is thus not possible to obtain A a based solely on the value of CTI. 

What is characteristic to the rhythmic oscillation in general is a steady zero-crossing 
pattern of the oscillation. Hence, a set of CTI's is considered of a origin if they correspond 
to continuous zero-crossing in the a rhythm range. Specifically, we consider a bigger set A 
using the following two criteria: (i) the CTI lies in the range of a rhythm range and (ii) the 
CTI's are from continuous zero-crossing. Let such CTI's be denoted by 



TOj > 1, i — 1, 2, • • •. Thus, Ai describes the i a event segment in the EEG record. The 
CTI's for the a dynamics is estimated by the union of A^s: 



Inevitably, some continuous zero-crossing could still come from the potential fractal com- 
ponent of the EEG. Hence, A a is only a subset of A. The "error" A\A a depends on the 
fractal property as well as the strength of the a rhythm. The complement C\A captures 
irregular zero-crossing that is characteristic to the fractal process 

For moderate a oscillation, EEG may contain other rhythmic components. The check of 
continuous zero-crossing must then be extended to a reasonably large scale range. To this 
end, we first consider the set of large CTI fluctuation: I\ = {r, > r u or Tj < r{\ where r u , T\ 




(la) 



A = uAi. 



(16) 
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are the upper and lower thresholds for defining the set of large CTI. Let /C denote the CTI's 
from continuous zero-crossings in (C\h)\A and let J = (C\/i)\(/C U.4.). Finally, the set of 
fractal CTI is obtained by 

T = h U Jo (lc) 

Note that this definition does not preclude fractal crossing in the a rhythm range. Note, for 
the a dominant case, that continuous crossing is mainly captured in A and, thus, /C ~ {0} 
and we recover T ~ C\A as shown above. If t u ^> /i(t) ^> 77, where /x(r) is the mean of 
{Tj}, any accidental fractal crossing included in /C U A will not introduce bias to the power 
law distribution of the fractal crossing (Fig. 2 below). 

B. Numerical Examples 

To demonstrate the above idea, we used Hilbert transform to construct the fractal time 
series given by the amplitude process of fBm, Ah(t). This is to mimic the fractal property 
reported in the a band-passed EEG 7 . Note that Ah(t) inherits the same scaling characteris- 
tics from B h (t). Hence, h = 2 — v holds, where v is estimated from the histogram of C\A. To 
define A, n ~ exp(— 5) and r u ~ exp(— 1) were used. These values are determined from the 
range of CTI of A h (t): r u ~ 0.8max(C), r\ ~ At, the sampling time. Fig. 2 shows the PDF 
estimate of C\A. It is seen that theoretical v values are verified before and after deleting A 
(Fig. 2). This should be the case since no band-limited component exists in A h (i). It thus 
establishes deleting A defined above does not affect the power law PDF of a pure fractal 
signal. 

To examine the influence from the band-limited rhythmic oscillation, Ah(t) in randomly 
selected time intervals of variable length were replaced by a narrow-band process x a (t) = 
M.{t)M{t) where M. = 1 + A h {t) models the fractal amplitude modulation 7 and the narrow- 
band M{t) is a sine wave of Gaussian amplitude X and frequency f ; i.e., X = N(l, o"x) and 
f = N(10,<Tf). Note that M(t) has a 10-Hz central frequency to mimic the a rhythm. To 
simulate the dominance of the band-limited oscillation, the probability of an interval being 
selected for x a (t) is four times of those for A h (t). In addition, the interval length for x a (t) 
is at least three times shorter than those for A h (t). The synthetic data so constructed is 
shown in Fig. 3a. 

In presenting our results, we keep the time unit in all figures so as to make easy reference 
to the narrow-band oscillation. Segments of CTI's before and after deleting A are shown in 
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Figs. 3b, 3c. The corresponding p(r) are estimated in Fig. 3d. It is clear that the narrow- 
band component x a can create significant bias in the otherwise power law PDF of Ah(t). 
The theoretical power law p(r) for A h (t) is correctly described after deleting A defined by 

(I)- 

III. EEG SCALING AND THE ALPHA DYNAMICS 
A. Scaling of the EEG Background Fluctuation 

We now apply the zero-crossing method to EEG records with varying degrees of a rhythm. 
These records were collected from six subjects (3 males, 3 females) of age 21 to 30 year- 
old (mean: ~24 who gave written consent to participate in the study. All subjects were 
instructed to maintain normal daily activity before participating in the 5-minute recording 
session. Surface scalp electrodes were attached according to the 10-20 international system 
at 01, 02, referencing to Cz. Two groups of data were taken: one in eyes open (EO) and 
one in eyes closed (EC). For EO, subjects were asked to direct their gaze at certain part 
of a shielded room to minimize eyes movements. For EC, no specific instruction was given 
to the subjects other than to relax and have their eyes closed. Output impedence from the 
recording system has been kept below 5kf2. The EEG was first band-passed from 0.1 to 70 
Hz and then digitized with a 12-bit A/D precision at 250 Hz (first four subjects) and 500 
Hz (last two subjects). 

In order to measure the strength of the a rhythm, we use the ratio of EEG signal power in 
the 8~12 Hz band to the full accessible frequency range, R a = J* g 12 S(f)df / J S(f)df (Fig. 4a). 
As expected, the R a is larger in EC than in EO due to the lack of visual stimulation in EC. 
Three of the six subjects (S2, S3, S4) are able to generate dominant a rhythm with large 
R a measure (> 0.45) in EC. It is important to note that low R a measure could mean very 
little a activity; e.g., the corresponding power spectra are given by power law with little 
identifiable feature in the a band. In contrast, there is always a very distinct "a peak" 
located at the 10Hz range for subjects showing large R a . 

The CTI PDF's of all EEG data sets are found to be of the power law form p(r) ~ r~ u . 
This indicates the fractal dynamics continues to exist in the brain state showing a rhythm. In 
particular, qualitatively different p(t)'s are found before and after deleting A from subjects 
showing dominant a rhythm (Fig. 4b), and almost the same p{r) from subjects showing 
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moderate to little a rhythm (Fig. 4c, also Fig. 2). For estimating the scaling in a dominant 
brain state, this means a more effective approach of using EEG zero-crossing than other 
amplitude-based methods such as the power spectrum and DFA (see also Fig. 1). 

Similar to the literature 8 , state dependence of the EEG fractal is observed: the v exponent 
is larger in EO than in EC. In addition, v and R a follows an inverse relation in both EC 
and EO (Fig. 5a). If one tentatively compares to the Hurst model using h — 2 — u, this 
result means a positive correlation between the EEG fractal background scaling and the 
underlying a rhythm. 

B. Alpha Interval and Organization of Complex Brain Dynamics 

The inverse relationship shown in Fig. 5a further suggests a relationship between the 
scaling property of the EEG fractal background and the a dynamics. 

Consider first the zero-crossing points of the EEG fractal background in the real line: 
If = {h, h, ■ ■ ■} where U locates the end points of the CTI defined in T . The set Lf can be 
obtained by intersecting the fractal portion of the EEG graph Xjr with the zero axis. Let 
dimi^ be the box-counting dimension for I jr. It is known from geometry that dim^ + l — 2 = 
dimi^ where the '1' on the left-hand-side of this equation is the dimension of the zero axis. 
By the Holder condition of the fractal function 19 , one can show dim^y < 2 — (2 — v) = v. 
Thus, dimi^ <v — \. 

Similarly, one can define the zero-crossing points of the a dynamics. In principle, this 
is obtained by intersecting the portion of the a oscillation in the EEG graph with the zero 
axis. In practice, we use the Aj's in (la) to approximate the a interval L; L = J2T=i t T k where 
the Li measures the "size" of the i th a event as it emerges in the foreground of the EEG. 
Unlike the in (la), which is narrowly distributed in the a rhythm range 18 , Lj can in 
general cover a much wider range of values. The end points locating the a intervals define 
the zero-crossing points I4 = {t±, t 2 , • ■ •}. 

Given the EEG fractal background established in the last section, we conjecture the same 
for the a interval and a power law PDF p(L) ~ . Denote the box-counting dimension 
of I4 by dim^; i.e., n(T) ~ T dimi A where n(T) is the number of zero-crossing points over 
the period T. The power law exponent <fi can be connected to dim^ via the relation 14 : 
n{T) = T — n(T) J T p(L)LdL where n(T) J T p(L)LdL on the right-hand-side estimates the 
average length of all a intervals over the period T. Substituting the power law p(L) ~ 
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yields dimi^ = 0—1. 

Since the point sets Lf and 1^ are embedded in the (ID) real line (time axis) and I^DI^ = 
0, one has 

= dirn^ + dim u - 1 < [v - 1) + (0 - 1) - 1. 

Hence, we have the inequality that connects the power law exponents of the EEG fractal 
background and the length of the a rhythm: 

z/ + 0>3. (3) 

Due to the little a rhythm in our EO data sets, the estimation of p(L) was suffered from 
poor statistics. For the EC data sets, both power law p(L) and (3) are verified (Fig. 5b). We 
will thus focus only on the EC data in this section. In Fig. 5c, an interesting transition is 
shown from the inequality (3) in a moderate cases (SI, S5, S6) to almost equality z/ + ~ 3 
in a dominant cases (S2, S3, S4). For z/ + ~ 3 (a dominant) the numerical values v ~ 1.25 
and ~ 1.75 are estimated. For v + > 3 (a moderate), larger v and are estimated. 

The case of ~ 1.75 in the a dominant EEG worths further discussion. In their model 
study of multi-trait evolution 14 , Boettcher and Paczuski (BP) predicted that evolution is a 
self-organized critical (SOC) process consisting of quiescent periods of all sizes, interspersed 
by short intervals of mutation event. These authors derived a similar power law scaling 
exponent of 1.75 for the distribution of the quiescent period and pointed out that their 
multi-trait model belongs to a new universality class. This universality condition prompted 
us to make the comparison of these two different phenomenologies. With the estimated 
~ 1.75, it is tempting to compare the a interval with the quiescent period in the BP 
model. This leads to the hypothesis of SOC dynamics in the "resting" state of the cortex. 

IV. DISCUSSION AND CONCLUSION 

The EEG background fluctuation coexisting with varying degrees of a rhythm is studied. 
It is important to note that our results directly address these two prominent features of the 
brain dynamics, rather than the EEG fluctuation in the a frequency band 7 ' 8 . Compared to 
other amplitude-based methods, we show that zero-crossing is more effective for studying 
the scaling in a dominant EEG. 

Our main result is the evidence and characterization of the coupling between the EEG 
background fluctuation and the a rhythm. Our findings can be summarized in two points: 
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(a) An inverse relationship between the EEG background scaling and the strength of a 
rhythm is observed, with a larger v exponent in EO for all subjects. Using the Hurst model 
tentatively (h — 2 — v), this implies a larger scaling exponent in EC compared to EO, or 
the trend towards more anti-persistent fluctuation in the a dominant brain state in EO. 

Similar state dependence propery has been reported in the past 3 ' 8 . Stam and de Bruin 8 
found similar result based on the correlation between the a band desynchronization and 
the decrease of the (DFA) scaling exponent in EO. An inverse relation was also reported by 
Moosmanns et al. 3 based on the blood oxygenation level dependence contrast as a measure 
of the brain metabolic activity. These authors concluded a dominant a rhythm associated 
with metabolic deactivation and desynchronization, a view may further be supported by the 
increase of local cortical activity 20 and information processing during EO 21 . 

(b) The inequality (3) characterizes the organization of the coupling. Since it is arrived 
based solely on the set-theoretic arguments, it is plausible that similar equations may exist 
between the EEG fractal background and other brain rhythms such as the 9, 5, (3, 7 waves. 

The observed transition from z/ + 0>3toz/ + 0~3ina dominant EEG, and the coinci- 
dence with the BP dynamics, imply a SOC state in the a dominant EEG. Hence, a strong a 
rhythm and the corresponding background fluctuation may represent two perspectives of the 
same dynamics. In general, this observation is in agreement with the suggestion of SOC of 
the EEG fractal dynamics 4 ' 7 ' 8 . However, our result differs in that we find indication of SOC 
only at the a dominant brain state, based mainly on the universality of the BP dynamics. 
For subjects showing little to moderate a rhythm, the unversality condition is no longer 
matched. We are not able to determine if a different universality class may be involved 
in these cases, nor can we ascertain a different theory for the observed fractal dynamics. 
Further studies on a larger population size and different physiological states are necessary 
to provide answers for these questions. They are the future work currently underway. 
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Figure Caption 



Fig. 1 (a) EEG records with moderate (top) and strong (bottom) a rhythm, and the 
corresponding (b) power spectral density functions and (c) the DFA result 11 F(l) ~ l b . The 
results for moderate (strong) a rhythm are the top (bottom) curves in (b) and (c). They 
correspond to, respectively, subjects S4 and S5 in eyes closed (Fig. 4). The solid line in (b) 
marks the frequency 10 Hz (log(10) ~ 2.3). The solid lines in (c) have the slopes ~1.21 and 
~0.51, respectively. 

Fig. 2 (a) An example of the CTI of A h (t), h = 0.3. (b) Log-log plot of p(r) for A . 3 (t) 
(top) and A , 8 (t) (bottom) before (open circles) and after (crosses) deleting A (see (1)). The 
axes are arbitrary. The solid lines are drawn with the theoretical slope — 1.2 (= 0.8 — 2) and 
— 1.7(= 0.3 — 2). The filled circles describe the zero-crossing PDF of a gaussian white noise, 
where no power law can be claimed 13 . 

Fig. 3 (a) A segment of the synthetic EEG. Two fractal periods are highlighted by horizontal 
bars in 161.5 ~ 162.5. (b) A segment of the set C. Note the concentration of t,- l ~ 0.05 
sec. (c) The set C\A where A is defined by (1). The horizontal lines mark the levels of 
t u ~ exp(— 2.5) and T\ ~ exp(— 5). (d) Log-log plot of p[f) before (connected dots) and 
after (heavy solid line) deleting A. The light solid line has the theoretic slope, -1.9, from the 
fractal component in synthetic data (A .i(t)). The vertical line marks the 10 Hz oscillation 18 
(r = l/20,log(r) ~ 2.9). 

Fig. 4 (a) The R a measure for the six subjects, (b) Log-log plot of p(r) for subject S4 in 
EC showing dominant a rhythm. A segment of the EEG for this subject has been shown 
in Fig. la. (c) Log-log plot of p{r) for a subject in EO showing little a rhythm (S6 in (a)). 
The close (open) circles show the p(r) before (after) deleting A and the solid lines are the 
regression lines (77 ~ exp(— 5) and t u ~ exp(— 1)). 

Fig. 5 (a) The inverse relationship between v and R a in EO (open circles) and EC (solid 
circles). The subject index (used in Fig. 4) is given next to the symbol, (b) Log-log plots 
of p{r) and p(L) for subjects showing dominant a rhythm in EC: S2 (in asterik and circle, 
resp.), S3 (in plus and triangle, resp.), S4 (in cross and square, resp.) Solid lines with slopes 
-1.25 and -1.75 are drawn, (c) The plot of v (solid circle), (open square) and z/ + (cross) 
for subjects in EC. R a in Fig. 4a is reproduced and referenced to the left y-axis. 
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